Duke preprint DU-TH-143 



Gluon Radiation and Coherent States in Ultrarelativistic Nuclear 

Collisions 



Sergei G. Matinyan^'^, Berndt Miiller^, and Dirk H. Rischke"'^ 

^ Department of Physics, Duke University 
Box 90305, Durham, North Carolina 27708-0305 

^ Yerevan Physics Institute, Yerevan 375036, Armenia 



We explore the correspondence between classical gluon radiation and quantum radiation 
in a coherent state for gluons produced in ultrarelativistic nuclear collisions. The expectation 
value of the invariant momentum distribution of gluons in the coherent state is found to agree 
with the gluon number distribution obtained classically from the solution of the Yang-Mills 
equations. A criterion for the applicability of the coherent state formalism to the problem of 
radiation in ultrarelativistic nucleus-nucleus collisions is discussed. This criterion is found to 
be fulfilled for midrapidity gluons with perturbative transverse momenta larger than ~ 1 — 2 
GeV and produced in collisions between valence partons. 
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I. INTRODUCTION 



The planned experiments at Brookhaven and CERN aimed at producing a dense quark-gluon plasma 
in highly energetic nucleus-nucleus collisions ||^ have generated interest in theoretical descriptions of 
the formation of such a state. One approach, the parton cascade model describes the equilibration 
process as a sequence of binary collisions of perturbative quanta in quantum chromodynamics (QCD). An 
unresolved issue of this model is how the decoherence of the initial nuclear parton distributions occurs. Two 
mechanisms have been proposed: the formation of a large number of mini-jets in the nuclear collision |^-^ 
and the radiation of soft gluons from the colliding valence quarks |^] . We are here concerned with the second 
mechanism. We note that it was recently argued in Q that the two approaches may be intimately related 
and lead to almost identical predictions in certain limits. 

The gluon radiation mechanism was originally |^] formulated in the framework of the classical approach 
proposed by McLerran and Venugopalan ||]. In their model, the small-a: gluon distribution of a large 
nucleus is described as arising from the classical gauge field created by a random distribution of color 
sources representing the valence quarks of the nucleons and moving along the light-cone, and by the quantum 
fiuctuations around this field. Kovchegov |0 has shown how the classical field can be obtained from a model 
of the fast moving nucleus as a collection of color dipoles and has identified the limitations of the classical 
approximation 

When two nuclei collide, a part of their classical fields is scattered on-shell, describing the emission of 
real gluons |7j. A detailed derivation of this process has recently been given to lowest and next-to- lowest 
order in the QCD coupling constant g [Q. Using semi-classical arguments, the radiated field energy can be 
related to the gluon multiplicity distribution as function of rapidity and transverse momentum. The result 
was shown to agree with the prediction of the perturbative, lowest order quantum calculation of soft gluon 
emission by colliding color charges JTst . 



In this work we demonstrate that classical gluon radiation as computed in ||7|,^ is closely related to 
quantum radiation of gluons in a coherent state, and in fact gives the same result for the gluon number 
distribution. In the coherent state, gluons are emitted in independent binary interactions between the 
colliding valence quarks of the two nuclei. Thus, if the number of emitted gluons is only a small fraction of 
the number of colliding color charges, correlations in the emission of individual gluons, such as arising from 
multiple collisions of the same charge, are negligible and the coherent state approach agrees with a quantum 
calculation of the multiple scattering process. If, on the other hand, a significant fraction of the colliding 
charges radiates, such correlations could become important and corrections to the coherent state calculation 
have to be accounted for. On the basis of this argument we estimate that for central Au + ylu-collisions 
at RHIC energies the coherent state approach is applicable to describe gluon radiation, if the gluons are 
produced in collisions between valence partons and at midrapidity with perturbative transverse momenta 
above about 1 — 2 GeV. 

The remainder of this paper is organized as follows. In Section II we briefly outline the results of Jl^ 
as far as they are of relevance to this work. Section III consists of two parts. In the first, the coherent 
state corresponding to gluon radiation from collisions of classical color charges is discussed. The second 
part interprets the results in terms of quantum diagrams. Conditions for the validity of the coherent state 
approach are discussed. Section IV contains the derivation of the event-averaged gluon number distribution 
in the coherent state formalism. The result is found to agree with that of Ref. We then generalize it to 
nuclear collisions at finite impact parameter. Section V concludes this work with a numerical evaluation of 
the gluon number distribution and a quantitative discussion of the validity of the coherent state approach 
for the description of gluon radiation in ultrarelativistic nuclear collisions. 

Our units are h ~ c = I, and the metric tensor is g'"' — diag(+, — , — , — ). Light-cone coordinates are 
defined as a± = (a°±a^)/\/2, dzp = d/dx±. The notation for 3-vectors is a = (a^,a^,a^) and for transverse 
vectors a = {a^,ay). 



II. CLASSICAL GLUON RADIATION IN ULTRARELATIVISTIC NUCLEAR COLLISIONS 

In this section, we briefiy summarize the results of fl^ , as far as they are required for the following. In 
covariant (Lorentz) gauge, d^A'^ = 0, the classical Yang-Mills equations, 

where = d^ — ig [A^ , • ] , F'"' = d^A'^ — d'^A'^ — ig [ , A'^], and J'^ is a classical source current, can 
be cast into the form 

□ ^J^'+ig[A^, d^A^' , (2) 

where □ is the d'Alembertian operator. In this form, the equations can be readily solved perturbatively 
order by order in the strong coupling constant, g. We expand 

oo oo 
fc=0 k=Q 

where AJ^"-*, J^"-* are the contributions of order g" to the gluon field and source current, respectively. [Eq. (^) 
takes into account that only odd integers of g occur in this expansion.] Then, the lowest and next-to-lowest 
order solutions obey 



□<) ^ J(i) ^ JW , (4a) 
□4^) = J(3) +^g[A('>, a. A« + Fi]} ] ^ J(3) . (4b) 



These equations are linear to each successive order in g and can therefore be solved with the method of 
Green functions: 

42fc+i) (^) ^ f ^4^/ ^^(^ _ j(2fc+i) , A: = 0, 1 . (5) 
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The retarded Green function reads in coordinate and momentum space [Q: 



Gr{x) = ^9{t)6{x^) , Gr{k) = - 



1 



(6) 



Let us now consider a collision of two nuclei with mass numbers Ai, A2, moving towards each other with 



ultrarelativistic velocities, wi 2 — ±1, along the z-axis. The nuclei are taken as ensembles of nucleons |10 
cf. Fig. |. 




FIG. 1. The nuclear collision as envisaged here. 



In order to simplify the color algebra we assume that each "nucleon" consists of a quark-antiquark pair. 
These valence quarks and antiquarks are confined inside the nucleons (visualized as spheres of equal radius 
in the rest frame of each nucleus). In order to construct the solution, nucleons inside the nucleus and valence 
charges inside the nucleons are assumed to be "frozen", i.e., they have definite light-cone (and transverse) 
coordinates. We assume the charges to move along recoilless trajectories, therefore, their coordinates will 
not change throughout the calculation. We label the coordinates of the quarks in nucleus 1 by X-i, x^, 
i = 1, . . . , Ai, and those of nucleus 2 by y+j, y j = 1, . . . , Antiquark coordinates follow this notation 
with an additional prime. 

Then, the lowest order classical current is a sum of the currents for each individual nucleus, as given in 
Eq. (3) of Re f. fl^ . The lowest order solution, A'i^\ and the associated field strength tensor, Pj^]) , were 
given in Ref. Jl^jTEqs. (4,5). This solution is identical to the one of the corresponding Abelian problem, i.e., 
the nuclei pass through each other without interacting. The fields generated by the valence charges simply 
superpose and no gluons are radiated. 

(3) 

With this lowest order solution and the classical current to next-to-lowest order in g, J/i , one can compute 
the next-to-lowest order solution, A^/^^ . The classical current jjj^"^ was obtained in from covariant current 
conservation, D^J^ — 0, and the assumption of recoilless trajectories, see Eq. (8) of The right-hand 



side of Eq. (4b) then reads in coordinate space 



where 

'^"(x- x,,y.^ = ^ [T^) (f;) (- 2^ ln(|x, - yJX) 5{x_ - x_,) e{x+ ~ y+,) S{x ~ x^ (8a) 

+ ln(|x - xJA) ln(|x -£^.|A)(9+(5(a;_ - a;_j)(5(a;+ - y+j)) , 
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P^'-ix; x,,y,) = ^ r"^ {Tt) (f/) ( 2^ \n{\x, - y^\X) e{x^ - x^,) 6{x+ - y+,) Six - yj 

- ln(|£-£jA) \n{\x-y \\)6{x- - x-i)d-5{x+ - y+j) 



—J ' 



f^^ix; x,,y,) = ^ r"^ {T^) (f/) - x^,) S{x+ - 



(8b) 



(8c) 



X ln(|x-y.|A) 



x — y. 

Mk-Xi\x) =^ 



Here, (Tf), (Tj) are color matrices which represent the color charge of the quarks in the color space of 
nucleon i of nucleus 1 and nucleon j of nucleus 2, while /"'"^ are the structure constants of SU{Nc). The 
antiquaries have the opposite color charge, — (T°), — (Tj"), which explains the relative signs in Eq. (|^), and 
ensures color neutrality of each nucleon. A is an infrared cut-off and acts as gauge parameter for the lowest 
order solution. 

The next-to-lowest order solution is then computed as 



-ik-: 



(27r)4 fc2 + ieko 



where 



and 



{k]Xi,yj) 



J-^°'{k]Xi,yj) 



J_ [k] Xi,yj) 



(27r) 
(27r)2 
(2n) 



(fq e 



-iir(x.-y,) 



{k^qf 



fc_ + ie q^ 
ik^ 



abcirpb\ /rf,c\ i{k+x-i+k-y+j-k-y ) 



d'qe 



-iq-{ 



qHk~qy 



(9) 
(10) 

(11a) 
(lib) 
(He) 



Note that an explicit expression for the solution Ajf''"(a;) was given in Rcf. |T^, Eqs. (21-24). The field (^) 
contains a piece associated with the change of color of a charge when it collides with the field of another 
charge. It arises from the pole k± — —ie in Eqs. (lla, llb| ) and thus does not correspond to radiated gluons, 
cf. also Eq. (25) in Jp2[ . The radiated gluon field, on the other hand, arises from the poles of the retarded 
propagator in Eq. (|ip^i.e., it corresponds to on-shell gluons, as one would expect. 



rad i^) — X/ ["^M rad i^l^iy Vj ) ~ rad (^'^^'i^ Vj ) ~ rad (^i , y'j ) + (x; X ■ , y'j ) 



where 



^!frad(^; ^^^V3) = ^(^ " ^^j" ) / ^k Uj'^^^^UJ, k; X,,y,) 



+ C.C. 



(12) 



(13) 



tij = [x^i + y+j) 1^/2 is the time when the collision between the quarks of nucleon i and j happens, and 
dk = (i^k/[(27r)^2aj], with lo = |k|. Obviously, this field vanishes prior to the collision. 



III. GLUON RADIATION AS A COHERENT STATE 
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A. Formalism 



In the preceding section and Refs. [[7yi2[|, the gluons produced in the colhsion of the two nuclei were treated 
as a classical field which solves the Yang-Mills equations. In contrast, in this section we want to describe 
the gluons as quanta radiated in the presence of a classical source, which we will take to be the current (Q). 
The way to solve a problem of this form can be found in standard textbooks jw) . 

Adapting the usual treatment to our case of a nuclear collision, we note that ai t —^ —oo the current (0) 
is zero and no real gluons are present. The gluonic Fock space "in" -state is therefore the physical vacuum, 
I ). When the two nuclei collide, the c-number source current produces real gluons such that at i ^ oo, 
the "out" -state of the system is the coherent state |p^ , |l5| 



\J) 



dkJ2Yl • e(k,A,a)c^(k,A,a) 



A=l a=l 



|0) , 



where 



(14) 



(15) 



is the on-shell Fourier transform of the classical current (0). It is identical with the current (^0|), taken with 
on-shell momenta fc^ — (lu, k), as implied by our notation. The e^(k. A, a) arc polarization vectors for (real) 
gluons with 3-momentum k, polarization A, and color a, obeying 



k ■ e(k, X,a) — , e(k. A, a) ■ e(k. A', a) ~ g 



XX' 



(16) 



The creation operators ct(k. A, a) and the corresponding annihilation operators c(k. A, a) for such gluons 
fulfill the commutation relation [Q: 

[c(k, A, a) , ct(k'. A', b)] = -2lu {2Trf g^^' Sab S{k - k') . (17) 

The first exponential factor in Eq. (^) ensures unitarity. It provides the normalization of the coherent state 
I J ) , with 



N 



dk 



j(3)"(w,k) • j'^'^>^{LU,-k) 



(3)a, 



(18) 



being the expectation value of the gluon number (with all possible colors and polarizations) in the coherent 
state. (The minus sign arises from our choice of metric, N is, of course, positive definite.) Since the color 
quantum number of final states in the collision will not be observed, it is not necessary to decompose the 
coherent state | J) in terms of states with good color quantum number jl^, i.e., for our purpose it is sufficient 
to consider the "conventional" coherent state (p^). 

The coherent state ( p^ is an eigenstate of the annihilation operator c(k. A, a) , 



c(k, X,a)\J) = -i J(3)'^(cj, k) • e(k. A, a) I J ) 



(19) 



Therefore, the expectation value of the number of gluons with momentum k, polarization A, and color a in 
the coherent state I J) is 



( J I ct(k. A, a) c(k. A, a) I J ) = J^^^^iu, k) • e(k. A, a) j(3)'^(w, k) • e(k. A, a) 
(no summation over a). Summing over colors and (transverse) polarizations, using 



^e^(k,A,a) e''(k,A,a) 



~9' 



{n ■ kf 



n ■ k 



(20) 



(21) 



with an arbitrary time- like unit vector n^, and current conservation 
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fc- j(3)'^(tj,k) =0 



(22) 



[this relation can be readily checked with the explicit form (|T^,|ll|) for the curr ent 4^^°], we obtain the 
invariant momentum distribution for the expectation value of the total gluon number in the coherent state 



dN 



dydP-k 2(27r) 
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— Y.Y.^J\^HK a, a) c(k, A, a) I J) ^ 



A=l a=l 



J(3)'^(a;,k)l • J(3)'^(cj,k) 



(23) 



Here y = ^ln[A:+/A:_] is the (longitudinal) rapidity of the gluons. Again, integrating over invariant momen- 
tum space dy d^k, we obtain the expectation value of the total gluon number N in the coherent state \ J), 
i.e., Eq. (H). 



B. Interpretation 

The coherent state | J), Eq. (^4|), is a superposition of 1-, 2-, . . . , n-, . . . gluon states. 



oo , ,„ n 



0) 



(24) 



dk ^ j(^)'''(wi,k,) • e(k,, Aj,a,)c'l'(k,, Aj,aj) 
(\0)^^ f dkY^T. ^~^'^"(^,k).e(k,A,a) |k,A,a) 

V A=la=l 

-^fd'kdk' J2 5] i^^^"(^,k) •e(k,A,a) J(3)^(w',k') •e(k',A',6) |k,A,a; k',A',6) + ... I . 

The amplitude for a final state of n gluons with momenta ki, . . . , k„, polarizations Ai, . . . , A„, and colors 
ai, . . . , a„ in the coherent state (e3) is 



A^(ki,Ai , ai; . . . ; k^, An, Q,^) = ( ki, Ai, ai; . . . ; k^, A^, Q,^ | 

■n. 

= e-^/2 n pj(3)-^(c^„k,).e(k„A„a,)' 



(25) 



i=l 



Thus, the probability to find n gluons in the final state regardless of spin, polarization, and color is 

Pn^ —f dki---dkn ^ ^ |Al(ki,Ai,ai;...; k„,A„,a„)|^ , (26) 



Ai,...,Ati ai,...,a„ 



where the prefactor takes into account that the gluons are indistinguishable Bose particles. A straightforward 
calculation yields the well-known result 



(27) 



where N is given by Eq. (|l|). The emission of gluons follows a Poisson probability distribution, or in 
other words, it happens in a statistically independent way. Although classically there are 4A1A2 collisions 
between the quarks and antiquarks, and all these collisions act as classical sources for the gluon field, quantum 
mechanically there is a certain probability that no gluon, one gluon, ... , n, ... gluons are emitted in the 
nuclear collision. This probability is given by Eq. (|27|). 
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p ^ p-q 
• f » > 



q-k 



p'-k+q 




E 




p'-k+q 



FIG. 2. Lowest order diagrams contributing to the amplitude for the emission of one gluon. 



The amplitude (|2^) for a state where n gluons are emitted during the collision can also be computed 
diagrammatically with the usual Feynman rules. Let us first focus on the case n — 1. To lowest order in 
5, the amplitude for the emission of one gluon with (on-shell) momentum A;^, polarization A, and color a 
is given by the diagrams in Fig. ^. In the recoilless limit for the quarks, a calculation of these diagrams 
with the usual Feynman rules yields the result — i J(^^°((jj, k) • e(k, A, a) [Note that in Ref. the gluon 

field to order g^, j4^^''°(a;), was calculated from the same set of diagrams. The only difference between 
that calculation and the present one is that now the emitted (on-shell) gluon contributes a factor (k, A, a) 
instead of the (retarded) gluon propagator times a phase e"**^'^ as in and there is no integration over 
the gluon 4-momentum k^.] 

To lowest order ing, the result for the perturbative amplitude for one-gluon emission is therefore identical 
to the amplitude ( P5| ) for n = 1. [To lowest order, the prefactor 1"^ ~ 1, since N ^ 0{g^).] This is 
another manifestation of the observation made in [^,|l2| that, to order g'^, classical and quantum calculation 
give the same result for gluon radiation. Note, however, that although classically there are AA1A2 collisions 
with gluon emission, the quantum diagram corresponds to only one single collision between any one charge 
in nucleus 1 and one in nucleus 2 and one single gluon emitted as a result of that collision. Of course, 

there are 4,AiA2 possible pairs of colliding charges, corresponding to the AA1A2 terms in J^^^"' , Eq. (p^). 
This increased probability for a collision is reflected in the gluon number distribution which also grows like 
4 cL Section IV. 

To higher order in g, there are the following type of quantum corrections to the one-gluon emission process: 

1. loops as shown in Fig. 

2. interactions of quark and antiquark inside a nucleon which participates in the gluon emission process, 
see for instance Fig. ||b, 

3. interactions between quarks and antiquarks of different nucleons in the same nucleus, cf. Fig. and 
finally 



4. interactions between quarks and antiquarks in different nuclei. Fig. 



abed 




FIG. 3. Higher order quantum corrections to the emission process. 
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The corrections (2), Fig. ^b, can be absorbed in the nucleon form factor (in the present approach, this 
form factor will appear after averaging over quark and antiquark coordinates inside a nucleon, cf. Section 
IV). Moreover, on the time scale of the nuclear collision gluon exchange inside a nucleon can be neglected. 
Similarly, we argue that corrections of type (3), Fig. |^, are nuclear structure effects and not important on 
the scales of interest here. More rigorously, it was shown in that diagrams of such type are exponentially 
suppressed in the covariant gauge. Corrections (1), Fig. and (4), Fig. are higher order quantum 
effects to the scattering process and suppressed by additional powers of g. Note that the corrections (2- 
4) induce correlations between the two charges directly involved in the gluon emission process and other 
charges. Following the above arguments, these types of correlations will be neglected. 

The amplitudes ( |25| ) for states containing n gluons are simply given by products of amplitudes for one-gluon 
states. Each one-gluon amplitude contains AA1A2 terms. A product of n such amplitudes, corresponding 
to emission of n gluons, contains (4 terms in total. Of these, (2Ai)! (2A2)! / {2Ai — n)\ {2A2 — n)\ 

terms correspond to processes where different charges collide once (in a given order) and emit a gluon, see 
e.g. Fig. ^a. These terms are identical to those one would obtain in a perturbative calculation in terms of 
diagrams, since these diagrams contain no multiple scatterings of a single charge. 



a b 




FIG. 4. 2-gluon emission processes, a) from collisions of different charges, b) from collisions where one charge 
collides twice. 

The remaining terms correspond to processes where a given charge collides and emits a gluon more than 
once, cf. Fig. In the coherent state, such a process is simply given by products of one-gluon amplitudes. 
Before and after each one-gluon emission process (as represented by the diagrams of Fig. ||), the charge 
is on-shell. On the other hand, in the corresponding diagrammatic calculation the charge is off-shell. In 
fact, since the nuclei are highly Lorentz-contracted and since the distance between subsequent collisions 
is inversely proportional to the corresponding Lorentz factor, the charge is rather strongly off-shell. The 
coherent state approach as discussed here fails to account for this and thus to give the correct description 
for the emission of more than one gluon from the same quark line. 

In other words, multiple collisions of the same charge induce correlations between the n emitted gluons. 
Such correlations are correctly accounted for in a quantum calculation in terms of diagrams. The coherent 
state approach, on the other hand, neglects these correlations and assumes that each individual gluon is 
emitted in an independent binary collision between color charges. This is what gives rise to the Poisson 
distribution (|27| ) for the emission of n gluons. Note that in addition to these correlations from multiple 
collisions, there are also higher order quantum corrections to the n-gluon emission process, which introduce 
correlations. These correspond to diagrams similar as in Fig. ||. According to the above discussion, these 
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correlations will be neglected here. 



(2A J ! (2A2) !/(2A, -n) !(2A2-n) ! 
(4A,A,)" 




In the following we want to derive a sufficient condition under which circumstances gluon emission from 
independent collisions of different charges dominates that from multiple collisions of the same charge. We 
argue that if this condition is fulfilled, the coherent state approach is applicable for the description of gluon 
radiation in ultrarelativistic nuclear collisions. We begin by noting that for small n and large Ai^ A2, the 
number of terms corresponding to multiple collisions of the same charge is negligible as compared to the 
number of terms where charges collide only once and emit a gluon. This is graphically shown in Fig. |^ where 
we plot the latter number (solid line) and the total number of terms (dashed). The right-hand panel shows 
the ratio of these numbers. Using Stirling's formula for large Ai,A2, this ratio can be approximated by 
exp[— n(n — 1) {Ai + A2)/4: A1A2], i.e., an approximate Gaussian behavior with n, as indeed seen in Fig. ^. 
Obviously, for n(ri — 1)^4 AiA2/{Ai + A2) the number of terms where charges collide only once and emit 
a gluon dominates the number of terms where gluons are emitted in multiple collisions of the same charge. 
In the situation of interest to us, A = Ai c:i A2 ^ 1, such that for the following discussion we may simplify 
this condition as <C 2 A. This condition is a sufficient criterion to neglect gluon emission from multiple 
collisions of the same charge. 

For 'n?<^2A, the n gluons are more likely to be produced in independent collisions of different charges 
than in multiple collisions of the same charge. In this case, we expect the description in terms of a coherent 
state to be (approximately) correct. For large ^ 2A, on the other hand, the emission of gluons from 
multiple collisions of the same charge is not negligible. The emission process for these gluons is incorrectly 
treated in the coherent state approach, and the discrepancy to the correct (quantum) result can in principle 
become large. Note, however, that a quantitative calculation is required to prove that this is indeed the case. 
It may be that even for large the deviation of the coherent state approach from the true answer remains 
small. Therefore, the criterion ^ 2 A is rather conservative. 

In order to minimize the deviation of the coherent state approach from the correct result, we have to require 
that the probability of occurrence of quantum states with a large number of gluons (which are incorrectly 
treated in the coherent state approach) is small. Since that probability follows the Poisson distribution 
( p7| ) in a coherent state, this happens exactly when the expectation value of the total gluon number in the 
coherent state, N, is small. Our criterion then becomes iV^ <C 2 A. Note that if ^ 2 A is fulfilled, 
it certainly holds for any subset of the total number N of radiated gluons as well. A subset of particular 
interest to us is, for instance, the number of gluons per rapidity, dN/dy. Moreover, even if it does not hold 
for the total gluon number, N, it may still be satisfied for a smaller subset such as dN/dy. This is obvious, 
since gluons in that subset are more likely to be less correlated. In Section V we shall discuss this criterion 
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in a more quantitative way for the situation of an ultrarelativistic nuclear collision. 

Note that this result is at first sight rather surprising. Naively, one would argue that the classical descrip- 
tion becomes better when the number of emitted gluons is large, such that the notion of a classical field 
becomes reasonable. As the description in terms of a coherent state is the closest analogue to the classical 
approach, one would expect a similar condition to hold in this case, too, for instance, that the expectation 
value of the number of gluons in the state | J) should be large. In contrast, from the preceding arguments 
one concludes that the coherent state approach is more viable for N"^ <^2 A. This apparent contradiction is 
resolved noting that not a large N is the essential criterion but rather that the emission of the n gluons in an 
n-gluon state happens in an uncorrelated fashion. This is exactly what gives rise to the Poisson distribution 
( p7| ) for the probability to find n gluons in the coherent state. If there are correlations in the emission of 
the gluons, the probability would no longer be a Poisson distribution, and the coherent state approach is 
no longer applicable. As discussed above, multiple collisions of the same charge are the major source of 
correlations and thus for potential deviations from the coherent state approach. The condition N'^ <^2A 
ensures that such correlations are negligible. 



IV. THE EVENT- AVERAGED GLUON NUMBER DISTRIBUTION 

The coherent state | J) depends on the initial color orientations of the quarks and antiquarks, (Tf ) , (Tj), 
and their coordinates Xi, x'^, yj, y'^ in a single collision event. These coordinates and color orientations are 
supposed to be fixed in the initial state. The coherent state is therefore a pure state. In other words, the 
phase information about the emitted gluons is complete, and all gluons are emitted coherently. This is 
certainly unphysical, because the emission of two gluons is incoherent if the spatial distance between their 
emission points is larger than their inverse momentum. In the classical approach of Refs. j^,^,^,^ as well 
as in the coherent state formulation discussed here, this incoherence is introduced a posteriori by averaging 
over all possible color orientations and coordinates. In the following, we therefore compute the respective 
average of the expectation value of the gluon number distribution (|2^) . We shall refer to this average as event 
average. In an extension of the treatment in we shall now also take a possible finite impact parameter 
h of the nuclear collision into account. 

Similar as in ||l2|, we introduce center-of-mass coordinates for nucleon i, (A_i, AJ, as measured in the 
center of mass of nucleus 1, and for nucleon j, (y+j, y^), as measured in the center of mass of nucleus 2, and 
relative coordinates Ax^i, Axj, Aj/_|_j, Ay . for the quarks and antiquarks inside each nucleon, as measured 
from the center of mass of the individual nucleon. Thus, if we assume that the center of mass of nucleus 1 
moves along the z-axis with velocity vi — 1, and the center of mass of nucleus 2 parallel to the z-axis at a 
transverse distance (impact parameter) b and with velocity V2 ~ —1, then 



Ax-^ , Ax-^ Ay+j , Ay+j 
, = — , y+j = r+j + , y^j = y+j — 



= + x^=X,-^, y^=h + Y,^ + ^, y^.=5 + y^.-^. (28b) 

In another generalization of the treatment in |T2[ we do no longer consider nucleons and nuclei to be 
"cylindrical" . The event average is therefore defined as 

d^XidX-^ fd'^{Axj2)d{Ax-^/2) 



X 



n 7i 47ri?3 / 3 V2 7 7i 47ra3 / 3 V2 7 
^ f d%dY,, r d\Ayj2)d[Ay,,/2) , ^ 



Here, R2 are the radii of nucleus 1 and 2, respectively, a is the nucleon radius, 7 the Lorentz factor in 
the center-of-mass frame of the nuclei, and the trace is taken over the color space of nucleon i and j. The 
integrations over "1" and "2" run over the volume of nucleus 1 and nucleus 2, respectively. Integrations over 
"i" and "j" run over the corresponding volumes of nucleon i and j. All volumes are Lorentz-contracted 
spheres. 
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We are now in the position to calculate the event- averaged expectation value for the gluon number dis- 
tribution, Eq. (p3|). If we follow the same steps as in |]l2| (and take "cylindrical" nucleons and nuclei) we 
indeed obtain the same answer, Eq. (39) of One of the assumptions of that derivation was, however, 
that the nuclei are (infinitely) large, cf. Eq. (35) of In that case, the notion of a finite impact parameter 
is obsolete. 

In the following we outline the derivation of the gluon number distribution in the case of a nuclear collision 
at finite impact parameter b and for spherical nucleons and nuclei. Inserting the current ( |l0|) into Eq. (|23|), 
the expectation value of the gluon number distribution in the coherent state | J ) becomes 



3" 



dN 

dy¥k ^ ^ 2(27r 

Ai A 

X 



abc rade 



d\^ d\. 



gl ■ 12^^ + fill ^ fll^' I2 ~ ll-' ll 



(30) 



^ g (rf ) (f;) {Ti) (fn V{k+,q^;x,)V{k^,k~ q^;y,)V*{k+,q^; Xk)V*ik^,k- q.^;yi) 



where k'^ is on-shell and 



V{k+,q\Xi) 



(31) 



We now perform the event average (p9|). The averaging over color is done utilizing tr [ {Tj') (TJ) (Tj}) (Tf) ] 
S^k S'"^6ji (5^74, and f'"' = N^N^ - 1), with the result 



d{N] 



dyd^k 2(27r)7 



li-l2^ +iiil-3.i^-i2-3.2^-ii 

ll ll (fc - {k-q2Y 



[V{k,q^,q^)), (32) 



where 



4=1 
A2 



g'(9i"£2'''^J^/^ _ e'''"'^^+^~ -■'^-i'^**-i"^-2^'^-j^^ + c c " (33) 
We now average over the light-cone variables Ax-i. The corresponding integration runs between if/"f = 

sin[2fc+//7] 2/ 



±V^a2- (Axi/2)2/y27 and gives 



//7 



//7 



Ax^, 



7 



(7 -> 00) 



(34) 



With an analogous relation for the average over Ay^j, the longitudinal momentum dependence in the 
phase factor drops out. The average over the light-cone variables X^i, Yj_j is even simpler, since V does not 
depend on them. The integration, however, yields a factor analogous to (|34). For the average over transverse 
coordinates we then use |1^ 



27ri?2/3 7|,|<« 



d^x e'^-^ \ 1-=- 



i?2 52^2 



M\i\R) 

\q\R 



cosi\q\R) 



^.h{\q\R) 

\q\R 



A{\q\R) , (35) 



where ji is a spherical Bessel function [note that for "cylindrical" nucleons and nuclei, A (a:) = 2Ji{x)/x 
with the result 



X A{\q^~q^\a)-A{\q^+q^\a)\ \ A{\q^ - q^\a) - A{\2k - q^ - q^\a) 



(36) 
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A(a;) is a rapidly decreasing function of its argument, thus the factors A(|g^ — q^\Ri), i = 1,2, limit the 
difference \q-^ — to a scale of order mm{TT / Ri , tt / R2} . That scale is much smaller than the typical scale 
of variation ^ 1/a of the remaining terms in (^ ) and in the integral over q^, q^ in Eq. (|3^). [An exception 

is the factor e*^-i~-2'' - which, for large impact parameters |&| ~ Ri, i = 1,2, then varies over the same small 
scale Tr/Ri, i = 1,2, and consequently must not be simply approximated by 1.] In these terms, we may 
therefore take to good approximation q^a q^. [This equality was achieved by quite similar means in [ p^ 



via the assumption that if i?i > i?2 > a, then A(|q^ - q \Ri) ^ Att/RI 5[q^ - q^), cf. Eq. (35) of |12[ .] 
The final result for the event-averaged gluon number distribution reads 

d{N) ^l^/T ,,,, 4 f F{\q\a) F{\k- q\a) 



dydPk (27r)6 N, ^ g2 J 1 ^2 (fc _ ^) 

with F{x) = 1 — A(2a;), and where we defined 



Ta,a, {b))^J^ e^L-'- A,A2 A{\l\R,) A(|l|i?2) . (38) 



While A in Eq. ( |37D (cf. the definition of F) stands for the nucleonic form factor [|13|, in ( p8[ ) it represents 
the nuclear form factor. 

We will now show that {TaiA2{1>) ) is the event- averaged number of binary nucleon-nucleon collisions per 
transverse area in an Ai + y4.2-collision at impact parameter b (the so-called nuclear thickness function). 
Note that for a central collision, uniform nuclear density distributions, and "cylindrical" nuclei, this factor is 
AiA2/ttR1 (for nucleus 1 being the larger of the two nuclei), and we thus re-obtain the answer Eqs. (38,39) 
of [|l^ [up to the modified form factor F{x)]. 

To prove this assertion, we revert the averaging over the center-of-mass coordinates of the nucleons, i.e., 

Ta.aM= fjikT.-'"-'-—-'^ (39) 



(2n) . . 



and introduce a factor of 



1 = J d^Xd^YS{X-Xi)SiY~Yj) (40) 

under the sum. This enables us to perform the /-integration. Introducing the baryon number distributions 
for nucleus 1 and 2, 

Ai A2 

PA,(x^i) = E'^(^-^»)^(^i-^^) ' PA.(i:,^2)-5]<5(i:-r,)'5(^2-^,) , (41) 

we realize that 

Ta^A2 ik) = Jd^xJ dZi PA, {X, Zi) j dZ2 PA, {X - b, Z2) (42) 

is the number of binary nucleon-nucleon collisions per transverse area for an Ai + yl2 -collision at impact 
parameter b. The event average of TAiA2{k) is simply 

( Ta,A2 (&) > - J d^X J dZi { PA, {X, Zi) > JdZ2{ PA2 {X - b, Z2) ) ■ (43) 

For Gaussian nuclear density distributions, pa(x) — A exp[— x^/i?^] (vri?^)^'^/^, we obtain the well-known 
result {Ta,A2{U)) = A1A2 exp[-feV(^? + ^1)] A(^? + ^1)- 

We have thus generalized the result of [T^ for the gluon number distribution to collisions at finite impact 
parameter. As expected, the gluon number decreases with the impact parameter according to the decrease 
of {TAiA2{]i)) with b. We also conclude from the results of this section that the classical solution |lj] and 
the coherent state approach yield the same expression for the gluon number distribution. 
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V. CONCLUSIONS 



In this work we have exhibited the relationship between classical gluon radiation in ultrarelativistic nu- 
clear collisions as calculated in JzlJ^JT^ and corresponding quantum radiation in a coherent state. We have 
demonstrated that the coherent state formalism yields the same result for the gluon number distribution 
as the classical approach. A condition for the coherent state approach to represent a valid description of 
gluon radiation in ultrarelativistic nuclear collisions was found to be that the expectation value of the gluon 
number in the coherent state has to be small compared to the number of sources in the (classical) source 
current, since then the emission of individual gluons happens in an uncorrelated fashion. The event-averaged 
version of this criterion reads (N)^ <C ZcS, where ZcS is the effective number of sources in the classical 
source current, which is, for the simple nuclear model considered here, equal to 2 A. 

This condition is similar to the criterion fs^ ^ ^ ^qcd ^'^^ applicability of the classical approach 
of Ref. since the latter states that the (transverse) area density of source charges should be large on 
transverse momentum scales of interest (such that the source can be described classically) while the number 
of field quanta ~ M^/fc^ associated with these sources must not be too large, i.e., it must be small as compared 
to //^/Aqqj-,. The problem is, however, that if one considers valence charges only, such as in the approach 
advertised here and in Rcfs. pO| , p^ , in the case of interest, i.e., for instance for a ylu -I- Ait-collision at RHIC 
energies, ^/s = 200 AGcV, ^^{N^ - 1)A/{2N^ iiR^) is only about (160 MeV)^, and thus of the order 
of Aqqq rather than much larger. It was therefore argued |^,|l8| that besides the valence charges considered 
here, one should also add quark and gluon charges arising from the nucleonic parton sea. 

On the other hand, in order to estimate whether the condition (iV)^ <C 2 A found here in the framework 
of the coherent state approach is quantitatively fulfilled, we compute the event-averaged gluon number 
distribution from Eq. (p7|). [We use the form factor for "cylindrical" nucleons, since in that case an analytical 
expression for the gluon number distribution is known [T^ ; the difference between the form factors in 
the cylindrical and spherical cases is small.] The result is shown in Fig. ^ for a central collision with 
Ai = A2 = 200. The strong coupling constant was taken as as = 3^/47r = 0.3, the nucleon radius (entering 
the nucleonic form factors) was assumed to be a = 0.8 fm. 
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1 2 3 4 5 

[GeV] 

FIG. 6. The gluon number distribution for a central A\ + A2-collision. 

As noted in [ p^ , the infrared behavior of the distribution is ~ 1/fc^, due to the nucleonic form factors, not 
^ log(fc^)/fc'' as claimed in [0J|]. Therefore, d{N)/dy diverges only logarithmically. This divergence is the 
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well-known infrared catastrophe [g^ . Physically, it will be regulated on a typical hadronic mass scale due to 
confinement. Here we simply introduce a lower transverse momentum cut-off to estimate d{N)/dy. 

Let us first choose the rather low \k^i^\ = Aqcd — 200 MeV. From a numerical integration of the spectrum 
in Fig. ^ up to |fc| = 10.2 GeV we obtain for the gluon rapidity density d{N)/dy ~ 140. The rapidity 
distribution is constant in y, since the radiating charges are assumed to move along the light-cone. However, 
since the classical approach fails anyway in the fragmentation regions, y ±t/cM p^ , one should not consider 
gluons other than such produced around midrapidity ?/ ~ 0. We shall therefore use the computed d{N)/dy 
(at midrapidity) to check the condition for the viability of the coherent state approach. As discussed in 
Section HI, it is always possible to use a subset for this purpose rather than the total number of gluons. We 
thus find that {d{N)/dyY ~ 20000 > 2A = 400, i.e., the criterion found in Section HI is not satisfied. In 
other words, the resulting gluon multiplicity is so large that a part of the gluons must have been produced 
in multiple collisions of the same charge. 

On the other hand, gluons with |fc| ~ Aqcd can certainly no longer be described in perturbative terms. If 
we consider an even smaller subset and take an infrared cut-off of 1 GeV, we obtain d{ N ) jdy ~ 26. For an 
even larger cut-off of 2 GeV ||], d{N)/dy ~ 8. For this cut-off, the gluon multiphcity is small enough that 
these gluons are likely to be produced in a statistically independent way. Thus, the coherent state approach 
to describe the production of such gluons seems justified. 

This result establishes that, although ^ ^qcD' i-^-' ™ principle violation of the fundamental assumption 
of the approach of Ref. 0, for practical purposes, i.e., for midrapidity gluons with perturbative transverse 
momenta |fc| > 1 — 2 GeV and produced in collisions between valence charges, a description in terms of 
a coherent state generated by a classical source of color charges is likely to be reasonable. The criterion 
/i^ ^ fc^ 3> Aq(^j-, is certainly a sufficient condition for the applicability of the classical approach. Here we 
have established that it might actually only be necessary to fulfill (7V)^ <C ^cff, with Z^s being the effective 
number of classical sources (and the total number of gluons squared, (N)'^, possibly replaced by a smaller 
subset of this number). If a quantum calculation of gluon emission from multiple collisions of the same 
valence charge proves that the deviation from the coherent state result for this process is small, then the 
range of applicability of the coherent state approach is, for practical purposes, even larger. Note that this is 
the case also for a smaller coupling constant. 

A gluon rapidity density oi d{N) jdy ^ 140 or even less for a central, ultrarelativistic collision oi A = 200- 
nuclei appears rather small. The reason is that, as already mentioned above, for A — 200 the average 
transverse color charge density squared = {N^ — 1)A/{2N^ ttR'^) is only about (160 MeV)^. If one 
adds the charge contained in the nucleonic parton sea |^,|l^, then fi ~ 400 MeV for collisions at RHIC 
energies |^ , which increases the gluon rapidity density by a factor of 40 (for fixed nuclear radius R) . For 
|fc^^.^| = 1 GeV, this leads to d{N)/dy ~ 1000, for = 2 GeV, we obtain d{N)/dy - 325. These values 

are well within in the established range |^J^ of several hundred gluons per unit rapidity, confirming the 
conclusion of Ref. 1^ that the classical description of gluon radiation and the conventional pQCD mini-jet 
approach give approximately the same results. 

One has to note, however, that increasing the effective charge of the classical current by a factor (3 (when 
adding the partonic sea) increases the gluon multiplicity by a factor (for fixed radii of the colliding nuclei) . 
Therefore, the criterion {d{N)/dy)'^ ^ ZcS, where now Z^s is the effective valence plus sea charge, may be 
violated in the process of adding charges from the nucleonic parton sea. In this case, multiple collisions of 
the same (valence or sea) charge become again important, and a more detailed quantum calculation of these 
processes is mandatory Q . 

Finally, we have generalized the result for the gluon number distribution found in | [T^ to nuclear collisions 
at finite impact parameter. As expected the number of binary nucleon-nucleon collisions per transverse 
area for given impact parameter, TAiA2ib), appears as a prefactor in the final result. 

In future applications it is important to incorporate non-Abelian effects on the dynamical evolution of the 
radiated gluons to study screening, damping, and, eventually, thermalization which leads to the formation 
of the proposed quark-gluon-plasma state of nuclear matter. 

ACKNOWLEDGMENTS 

D.H.R. thanks Miklos Gyulassy and Yuri V. Kovchegov for valuable discussions. The authors are indebted 
to Yuri V. Kovchegov for a critical reading of the manuscript. This work was supported in part by the U.S. 



14 



Department of Energy under Contract No. DE-FG02-96ER-40945. 



[1] 
[2] 
[3] 

[4] 



[5] 

[6] 

[7] 
[8] 
[9] 
[10] 

[11] 
[12] 
[13] 
[14] 
[15] 
[16] 
[17] 

[18] 

[19] 



B. Miiller, Rep. Prog. Phys. 58, 611 (1995). 

J.W. Harris and B. Miiller, Annu. Rev. Nucl. Part. Sci. 46, 71 (1996). 
K. Geiger and B. Miiller, Nucl. Phys. B369, 600 (1992); 
K. Geiger, Phys. Rep. 258, 237 (1995). 

K. Kajantie, P.V. Landshoff, and J. Lindfors, Phys. Rev. Lett. 59, 2527 (1987); 

J.P. Blaizot and A.H. Mueller, Nucl. Phys. B289, 847 (1987); 

K.J. Eskola, K. Kajantie, and J. Lindfors, Nucl. Phys. B323, 37 (1989). 

X.N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 (1991); Phys. Rev. D 45, 844 (1992); Phys. Rev. Lett. 68, 

1480 (1992); Phys. Lett. B 282, 466 (1992); Comp. Phys. Comm. 83, 307 (1994). 

K.J. Eskola and M. Gyulassy, Phys. Rev. C 47, 2329 (1993); 

K.J. Eskola, A.V. Leonidov, and P.V. Ruuskanen, Report No. hep-ph/9606406 . 

A. Kovner, L. McLerran, and H. Weigert, Phys. Rev. D 52, 6231 (1995); 52, 3809 (1995). 

M. Gyulassy and L. McLerran, Report No. tiucl-th/9704034 , 

L. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994); 49, 3352 (1994); 50, 2225 (1994). 
Yu.V. Kovchegov, Phys. Rev. D 54, 5463 (1996). 
Yu.V. Kovchegov, Phys. Rev. D 55, 5445 (1997).^ 



Yu.V. Kovchegov and D.H. Rischke, Report No. |hep-ph/970420l| (submitted to Phys. Rev. C). 
J.F. Gunion and G. Bertsch, Phys. Rev. D 25, 746 (1982). 

C. Itzykson and J.-B. Zuber, Quantum Field Theory, (McGraw-Hill, New York, 1985). 
C. Slotta and U. Heinz, Phys. Lett. B 391, 469 (1997). 

K.-E. Eriksson, N. Mukunda, and B.-S. Skagerstam, Phys. Rev. D 24, 2615 (1981). 

LS. Gradshteyn and LM. Ryzhik, Table of Integrals, Series, and Products , (Academic Press, San Diego, 1980) 
Eqs. (6.567.1) and (8.464.3). 

J. Jalilian-Marian, A. Kovner, L. McLerran, an d H. Weigert, Report No. hep-ph/9606337; 

J. Jalilian-Marian, Report No. hep-ph/9609350. 

S.G. Matinyan, B. Miiller, and D.H. Rischke (in preparation). 



15 



